ggplot(D)+
stat_smooth(aes(x=Shade.Index, y=Fv.Fm), formula = y ~ x, method = "lm", se = TRUE, color="grey", linetype= "solid", size=0.5)+
geom_point(aes(x=Shade.Index, y=Fv.Fm, fill=Veg.Type, size =PDIR), shape= 21, color = "black") +
scale_size(name = "PDIR", range = c(0.25, 6)) +
stat_smooth(aes(x=Shade.Index, y=Fv.Fm, color=Veg.Type),
formula = y ~ x, method = "lm", se = FALSE)+
labs(y='\u2190 Stress (Fv/Fm)', x="Microscale shade time (% of the year)")+
scale_y_continuous(breaks = seq(0.1, 0.75, 0.1), limits=c(0.1,0.75))+
#annotate(geom="text",label= "R-sq = 0.082\nP < 0.0051", x=5, y=0.7)+
#scale_color_manual(values =c("#E6AB02", "firebrick3", "dodgerblue3"))+
scale_fill_manual(values =c("#E6AB02", "firebrick3", "dodgerblue3"), name= "Elevation\n(life zone site)", guide = guide_legend(reverse=TRUE))+
scale_color_manual(values =c("#E6AB02", "firebrick3", "dodgerblue3"), name= "Life zone", guide="none")+
#guides(color = "none")#guide="none")+
theme_linedraw()+
theme(axis.text = element_text(size=12), axis.title = element_text(size=15)) #text = element_text(size=15))#must place after chart theme